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Abstract 

We consider a spin ladder model which is known to have matrix product states as exact 
ground states with spin liquid characteristics. The model has two critical-point transitions 
at the parameter values u = and oo. We study the variation of entanglement and fidelity 
measures in the ground states as a function of u and specially look for signatures of quantum 
phase transitions at u = and oo. The two different entanglement measures used are S(i) (the 
single-site von Neumann entropy) and S(i,j) (the two-body entanglement). At the quantum 
critical point (QCP) u = oo, the entanglement measure E [= S(i), S(i,j)] vanishes but remains 
non-zero at the other QCP u = 0. The first and second derivative of E with respect to the 
parameter u and the entanglement length associated with S(i,j) are further calculated to 
identify special features, if any, near the QCPs. We further determine the GS fidelity F and a 
quantity ln\D\ related to the second derivative of F and show that these quantities calculated 
for finite-sized systems are good indicators of QPTs occurring in the infinite system. 

I. INTRODUCTION 

In recent years, quantum phase transitions (QPTs) in many-particle systems have been extensively 
investigated using well-known quantum information theoretic measures. QPTs which are solely 
driven by quantum fluctuations occur at zero temperature when some parameter, either external or 
intrinsic to the Hamiltonian, is tuned to a special value termed the transition point [l]. In the case 
of second-order QPTs (critical-point transitions) , a diverging length scale governs the physics near 
a quantum critical point (QCP). Usually, the correlation length associated with specific correlation 
functions diverges as the QCP is approached and the ground state properties develop non-analytic 
features. In this context, it is pertinent to ask whether the quantum correlations associated with 
entanglement are good indicators of QPTs. A number of entanglement measures have so far been 
identified which show special features close to the transition points of QPTs occuring in spin systems 
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El [H [5J [6l \7\. It has been shown that, in general, a first order QPT, linked to a discontinuity in the 
first derivative of the ground state energy, is signalled by a discontinuity in a bipartite entanglement 
measure such as negativity and concurrence [H [H [10] and a discontinuity or a divergence in the first 
derivative of the same entanglement measure marks a second order phase transition characterized 
by a discontinuity or a divergence in the second derivative of the ground state energy |3j. The 
entropy of entanglement of a block of L contiguous spins in a chain with the rest of the system has 
been shown to diverge logarithmically with L near the QCP [4]. 

The typical length scale over which a particular entanglement measure decays defines the entan- 
glement length (EL) . A number of entanglement measures characterized by a diverging EL close to 
a QCP have been proposed to date [Hill]. One of these, the two-body entanglement S(i,j) which 
estimates the amount of non-local correlations between a pair of separated spins at sites i and j 
and the rest of the spins, is given by the von Neumann entropy 

= -Trp(i,j)log 2 p(i,j) (1) 

where p(i,j) is the two-site reduced density matrix obtained by tracing out the spins except the 
ones at sites i and j from the full density matrix. When the system is translationally invariant, S 
depends only on the separation n =| j — i | and can be expressed in terms of the spin correlation 
functions in the large n limit. Away from the critical point, S(i,j) saturates over a length scale £e, 
the EL, as n increases. Close to the QCP and for large n, we have 

S(n)-S{oo) ~ A{n)e~^ (2) 

S(i,j) captures the long-range correlations at a QCP if A(n) has a power-law decay as a function 
of n along with a divergent This is true for spin models such as the S = 5 exactly solvable 
anisotropic XY model in a transverse magnetic field [7] . The EL is found to diverge with the same 
critical exponent as the correlation length at the QCP. S(i,j) and its first derivative have been 
found to develop special features in the vicinity of the QCP [HI [HI [20]. The single-site von 
Neumann entropy (a measure of the entanglement of a single spin with the rest of the system) 

S(i) = -Trp(i)log 2 p(i) (3) 

is also known to be a good indicator of QPTs El [TBI [19]. In Eq. (3), p(i) is the single-site reduced 
density matrix. 

The exploration of the entanglement properties of the ground state of a number of spin— 1 
Hamiltonians (the AKLT model is an example) has been carried out using both analytical and 
numerical techniques [HI 22 02] . Certain spin-1 and generalized spin- 5 ladder model systems are 
known to have matrix product (MP) states as exact ground states [Ml [HI [16] . The MP states 
are finitely correlated states with short-ranged spin-spin correlations, may have hidden topological 
order and have gapped excitation spectra. The second order transitions in these so-called finitely 
correlated MP states belong to the class of generalized QPTs (the definition encompasses the tran- 
sitions marked by a non-analyticity in any observable of the system) [TB] which differ from the 
conventional QPTs in some important aspects. The spin correlation function in both the cases is of 
the form Ac e~^c for large n. In the case of MP states, Ac vanishes at the transition point though 
the correlation length £c blows up as the transition point is approached. In the case of a con- 
ventional QCP, the correlation function has a power-law decay close to the QCP. A distinguishing 
feature of QPTs in MP states arises from the fact that the ground state energy density is analytic 
for all values of the control parameter. A critical point transition is, however, still signalled by a 
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diverging correlation length and the vanishing of an energy gap. The MP states have been used as 
trial wave functions for a number of standard spin models and provide the basis for the well-known 
density matrix renormalization group (DMRG) method leading to several interesting developments 
in quantum information theory [17]. The MP states further serve as candidate systems for the 
study of unconventional QPTs. 

Recently, ground state fidelity has been proposed to provide a signature of QPTs [2Tj [22] and 
the usefulness of the measure has been explored in a number of studies [23, 24, 25, 26, 27, 28, 29J. 
Fidelity, a concept borrowed from quantum information theory, is defined as the overlap modulus 
between ground states corresponding to slightly different Hamiltonian parameters. The advantage 
of using this measure is that it characterizes QPTs without needing any a priori knowledge of 
the order parameter and the symmetries of the system. The fidelity typically drops in an abrupt 
manner at a transition point indicating a dramatic change in the nature of the ground state wave 
function. A QCP is characterized by the vanishing of the single particle excitation gap. In Ref . [27J, 
an explicit connection between the vanishing of the gap and the fidelity drop has been established. 
Cozzini et al. [23j tested the validity of the fidelity approach for probing QPTs in MP states and also 
studied the finite size scaling of the fidelity derivative establishing its relevance in extracting critical 
exponents. The QPT in the Bose-Hubbard model which is difficult to detect using conventional 
entanglement measures has been correctly predicted using the fidelity measure [25]. Chen et al. [29] 
have shown that the fidelity of the first excited state and not the ground state, is the appropriate 
quantity to signal QPTs in models such as the antiferromagnetic (AFM) Heisenberg spin chain with 
nearest-neighbour as well as next-nearest-neighbour interactions. 

In this paper, we study a S = \ ladder model with MP states as exact ground states [30J. 
The model has an interesting phase diagram with two critical point transitions. We explore the 
properties of the ground state using two different entanglement measures, namely, the single-site 
and the two-body entanglement. The major motivation is to identify distinctive features, if any, in 
the entanglement measures close to the QCPs. We look at the same QPTs in the light of fidelity 
analysis and show that the fidelity T of the ground state is an efficient indicator of the QPTs. The 
quantity In \D(u)\, related to the second derivative of T , also yields useful information regarding 
the QPTs. We apply the idea of average entanglement [31] to take care of the two-fold degeneracy 
of the ground state of the model. 

II. ENTANGLEMENT MEASURES 

We consider a general ladder model proposed by Kolezhuk et al. [30] which is described by a 
Hamiltonian of the general form 

H = 'y^[J(S± t jSl J j+l + S2,jS2,j+l) + JrSljSsj + V(Si,jSij + i)(S2,jS2,j + l) 

+ Jd(Si,jS2,j+i + S^.jSij+i) + X{(5i j S , 2 J+ i)(S , 2jS'i j+ i) - (SijS2,j)(Sij+iS2,j+i)] (4) 

where the indices 1 and 2 distinguish the lower and upper legs of the ladder and i labels the rungs. 
The ground state ji/'o (u, u)) has the following MP form 

|V>o (u, it)) = -j==Tr {0i(u).02(m) 92N-i(u).g2N{u)) (5) 
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where 

and ATc (= (uu + 3) 2N + 3 (im — l) 2N ) is the normalization factor. Here \s) i is the singlet state 
and \tfj_) with /x = +1, and —1 are the triplet states of the i-th rung with S z = + 1, and — 1, 
respectively. 2N is the total number of rungs (with periodic boundary conditions) and u, u are free 
parameters. For u ^ u, \4>o (u, u)} is dimerized and doubly degenerate as the translation of the 
rungs by one unit leads to a different state with the same energy. 

It is convenient to rewrite the Hamiltonian (5) as a sum of identical local terms which couple 
only neighbouring rungs, H — ^(^M+i — Eq). The value of E Q is adjusted to make |V>o) a zero- 
energy ground state which requires the following conditions to be satisfied, (i) All elements of the 
two matrix products <?i(u).<?i+i(u) and gi(u).gi+i(u) have to be zero-energy eigenstates of hij+i. 
(ii) The other eigenstates of should have positive energy. The two conditions are satisfied 

when hi t i + i has the structure 

./ 

h i;i+ i= ^2 ^2 e j\i ) JM) (iPjm\ (7) 

,7=0,1,2 M=-,J 

where the eigenvalues ej > and |^jm)'s are the components of the positive-energy multiplets 
constructed from the states of the two-rung plaquette (i, i + 1): 

\tpoo) = [3 + {uu) 2 ]-i{V3\ss) +uu\tt) J=0 } 

Hi) = [2 + (m-H\st) + \ts) + f\tt) J=1 } (8) 

u-\-u 
V2 



The notation \tt) J=1 has been used to describe states with the total spin J = 1 constructed from 
two triplets on rungs i and i + 1, etc. We obtain the connections between the parameters J, J r , 
Jd, V and K of Eq. (4), the local eigenvalues ej and the singlet weight parameters u, u of the 
ground state wave function by claiming that the structure (7) is compatible with Eq. (4). The 
model we study in this paper is a special case of the three types of solutions obtained from the 
above-mentioned relationships. In this case, Jd — 0, K ^ and 

u = -u, K = J r = eo (u2 '^ u2+3 \ J d = 0, 

y = eo (6u*+2u'+9) i j = (ttVlOn'+S) | (9) 

_ (3u 4 + 14« 2 + 15) (5m 4 + 18« 2 +9) 
£ 1 — e g , £2 — £0 § 

As pointed out in Ref. [30J , the one-parameter model undergoes two second-order phase tran- 
sitions, one at u = and the other at u = oo. At u = 0, the ground state undergoes a transition 
from the dimerized phase to the Haldane phase. The effective Hamiltonian describing this phase 
is that of the S = 1 AKLT chain. At u = oo, the transition is to a phase in which the ground 
state is a product of singlet bonds on the rungs. The transitions at u — and oo are marked by 
the vanishing of the singlet and triplet gaps, respectively, in the excitation spectrum [30] • The 
ground state is spontaneously dimerized everywhere except at the critical points. In the MP 
formalism, it is straightforward to calculate the spin-spin and dimer-dimer correlation functions 
C s (n) = (Sfj i+n )and C D {n) = (A A+n) where A = S M .(S M+ i - S M _i). The dimer corre- 
lations are long-ranged and vanish as u — > 0, oo but with no exponential tail. The spin correlation 
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length is finite at the AKLT point u — 0, becomes zero at u — 1 and diverges as u — > oo. There 
is, however, no development of long-range spin order since the amplitude of the spin correlations 
becomes zero in this limit. The doubly-degenerate spontaneously dimerized phase which prevails 
away from the critical points exhibits non-Haldane spin liquid properties. The elementary excita- 
tion is of a novel type, a pair of propagating triplet or singlet solitons connecting two spontaneously 
dimerized ground states [30]. In the Haldane phase, the elementary excitation has the character of 
a magnon. 

Using the transfer matrix (TM) method, we now study the entanglement properties of the MP 
ground state [Eq. (9)]. The state is two-fold degenerate as the ground-state energy per rung 
Eq = A (7u 4 + 22u 2 + 19) does not depend on the sign of u. The two ground states obtained 
from Eq. (6) 



l^ 1 ) = /m r , Tr {9l(.u)-92(-u) 92N-l{u).92N{-u)} 

\ip2) = 7= Tr{ gi {~u).g 2 (u) 92N-i(-u).g 2 N(u)} 



(10) 



are asymptotically orthogonal in the thermodynamic limit (TDL) N — > oo, i.e., the overlap 
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(V'llV^) = ( i "+ 3)Vjv + ^ [ ~5 _ 1 1 a n < 1 for finite N and vanishes in the limit TV — » oo. No(u) 
(u 2 + 3) 2N + 3 (u 2 — l) 2N ] is the normalization factor. We construct a pair of orthogonal degenerate 
ground states applying the usual Gram-Schmidt procedure 

\<h) = \1>i) 

\cp 2 )= i (|^ 2 )-(Vi|^)|Vi)) ( u ) 

VJf 



with N = 1 — KV'ilV'a) | 2 - An arbitrary superposition of the two degenerate ground states is also a 
valid ground state. We apply the idea of average entanglement [31], i.e., calculate the entanglement 
content of a general state (an arbitrary superposition of basis states) and then calculate its average 
value over the whole of parameter space (the coefficients of the basis-state expansion constitute the 
parameters) 

f dn(p 1 ,p2,...)\E(p 1 ,p 2 ,...)\ 

J dfl[puP2, ■■■) 

where J d/j,(pi,p2, •••) is the Haar measure associated with the parametrization p%,p2,..., which 
is invariant under unitary operations. The normalization condition restricts the values of the 
coefficients so that the parameter space is associated with a compact hyper-surface. In the case of 
a double degenerate ground state, a general state is a superposition of two states 

\4>„) = a ](/)!) +b\(f> 2 ) (13) 

with the restriction \a\ 2 + \b\ 2 = 1. The corresponding parameter space is a 3 — D sphere S 3 . The 
one-rung reduced density matrix p(i) (Eq. (3)) is obtained by tracing out all the rungs except the 
i-th one from the ground state density matrix p — \<j> s )((f> s \ . From Eq. (13) 

P(i) - Trl, L \<p s ) = Tr\^ L {\a\ 2 |0 X > + \b\ 2 \cf> 2 ) (0 2 | + ab* |^) (0 2 | + a*b (14) 

With the help of standard TM calculations [19] , one obtains a form for p{i) which is found to be 
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o^oo 

^T3 



(15) 



independent of the parameters a and b in the TDL 

/ 

Pii) = 

V ^/ 
in the |i±i,o, s) basis. The single-rung entanglement is obtained as 

S ( l ) = -Y—Ti [(" 2 + 3 )^°32 (" 2 + 3) - u 2 log 2 u 2 ] (16) 

Entanglement average, as defined in Eq. (12), is required for finite-sized systems. In the TDL, 
such averaging is not necessary as p(i) [Eq. (14)] is independent of a and b (\a\ 2 + \b\ 2 = 1). The 
variations of S(i) and its first derivative with respect to u have been shown in Fig. 1 (top) and 
(bottom) respectively. S(i) has the value logi 3 at the critical point u = (the AKLT point) as 
expected, increases as u is increased from zero before it reaches its maximum possible value of 2 at 
u = 1. Then it decreases with increasing u and vanishes at the other transition point u — oo (Fig. 
2). In the rung-singlet phase, each pair of spins in a rung forms a singlet to become maximally 
entangled with each other and completely unentangled with the rest of the system. The plots are 
expectedly symmetric about the point u = 0. 

The two-rung reduced density matrix p(i,j) can be calculated in the same manner. p(i,j) is 
given by 

p(i,j)=Tr\l L \<t> s ){<j> s \ (17) 

where the trace is taken over all the rungs except the i-th and j-th ones. From the usual TM 
calculations , we obtain p(i,j), in the TDL, as a 16 x 16 matrix in block-diagonal form. From (1) 
and (17), the two-body entanglement is 

s(hj) = -J2 x * lo ^ x * ( 18 ) 

i 

Aj's being the eigenvalues of p(i,j). Figure 3 shows the variation of the average S(i,j) (top) and 
its first derivative (bottom) with u for n = 1000. S(i,j) behaves in a similar manner as S(i). It has 
the value 2 log2% at the QCP u — 0, it then increases with u to attain the peak value 4 at u = 1 and 
when u is increased further, S(i,j) decreases and falls to zero (Fig. 4) as we approach the QCP 
u = oo . The first derivatives of S(i) and S(i,j), instead of showing any non-analyticity, fall sharply 
to zero at both the QCPs. The first derivatives are also zero at u = 1 where the entanglement 
measures have the maximum value. The second derivatives of S(i) and S(i,j) are logarithmically 
divergent at both the QCPs u — and oo (as can be seen in the insets of Figs. 1 — 4). Both the 
measures S(i) and S(i,j) vanish at u = oo and are non-zero elsewhere, they thus behave as an 
order parameter decreasing to zero value at the QCP u = oo with transition to the rung-singlet 
phase. The measures, however, do not have the character of an order parameter for the transition 
at u = from the dimerized to the Haldane phase. 

We next calculate the EL, £e, associated with the entanglement measure S(i,j) . Close to either 
of the QCPs and in the limit of large n, we have S(n = \j — i\) — S(oo) <~ A e e . The pre-factor 
A e is found to remain finite and non-zero at the transition point u = but it vanishes at u = oo . 
The EL £e is given by 

& - iski (19) 
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We rewrite £e as a function of -7, i.e., £e = rnr, 1 ^2 and study its behaviour near u — 00 i.e., 

~ = 0. Fig. 5 shows the variation of £e with respect to u and ~. We find that £g is finite at the 

critical point u = but it diverges as u — > oo with the critical exponent ^ = 2 as £e ~ (i) for 
~ ~ 0. The spin-spin correlation function Cs{n) —< Sf i Sf i+n > can be calculated in the TM 
formalism as |30j 

C s (n) = (u 2 + 3) _1 (z + z_) n (<5 n , 2 fe - z- <5 n ,2ft+i) , on x 
z± = ( U ±l) 2 /(7i 2 + 3) ( ^ Uj 



Close to the QCP u — oo, ~ £c/2 so that both £e and £c diverge with the same exponent 
v = 2. 

II. GROUND STATE FIDELITY JF(w, <J) 

We now investigate the behaviour of fidelity near the same pair of QCPs. The average fidelity, in 
analogy to (12), is 

T _ J d(J,(p 1 ,p 2 ,...)J r (pi,p 2 ,-) , s 

J~ av — , , \ ^J-J 

OM(Pl>P2, ■■■) 

The overlap between two general ground states, |0(ui)) and |0(u2)) (see Eq.(13)), corresponding 
to two different values of the control parameter is given by 

T( Ul ,u 2 ) = (0(u x )| 0(u 2 )> = |a| 2 (0i| 0i> + H 2 (02 1 02) + a6* (0i| 2 ) + a*& (02 1 0i) (22) 
T{u\, u 2 ) (averaged over the {a, b} ) can be expressed in terms of the eigenvalues of the TM [23] as 

T( Ul ,U 2 ) = -i 1 [(1 + l+P("i)P("2) u ( + 3 )2N + 3 

(^2-i) 2Jv >- A / ( y; ( "if ,, {(^2-3)^+3(^+1)^}] 1 J 

V( 1_ P 2 («l)P 2 («2)) 

3 (u 2 + l) 2iY -t-(u 2 — 3) 2JV 

where p(u) = ^+3)'^ + 3(^-1)'^ ' 6 (top) shows the nature of the variation of u + S), 

(overlap of the states corresponding to two close points in the control parameter space separated 
by a small variation) with u and N in the neighbourhood of the critical point u — for 8 = 
.001. A straightforward calculation reveals that for large values of N and for m^O, T{u\,u 2 ) ~ 

(a(ui, u 2 )) N , where a(u\, u 2 ) — "^l^f" 1 ," 2 ^^ ■ a(u, u + 5) < 1 and it has a sharp dip at u = 0. 
Thus away from the critical point, u + S) decreases exponentially with N and vanishes in the 
TDL for any fixed value of u and S, but we observe from Fig. 6 (top) that J-(u, u + S) decreases at 
a much enhanced rate when the QCP is approached. Intuitively, the rate of orthogonality, i.e., the 
rate at which the "distance" between the ground states corresponding to two neighbouring points 
of the parameter space becomes maximal, should diverge in the proximity of a QPT. It is thus 
sensible to relate the degree of criticality to the derivative of the fidelity function. Cozzini et al [23] 
have proposed a general expression for the quantity relevant in this case 

D(u) = -d Ul d u JnF(ui,u 2 ) | Ul =„ 2 =« (24) 



where F(u\,u 2 ) = J N§{u\)Nq(u 2 )T(u\, u 2 ). In the large N limit and for 11 ^ 0, one can easily 

N 



check that D(u) ~ CTXgp ■ Thus in the plots (Fig. 6 (bottom) ) showing the variation of ln\D(u)\ 
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with u for different values of N, we observe that the rate at which ln\D(u)\ increases with u is 
heightened in the proximity of the QCP u = 0. To repeat the whole analysis for the other critical 
point u = oo, we express F(ui,u 2 ) as a function of U\ = — and u 2 = —, For very large N, 

T'(ui,u 2 ) ~ (a'(ui, u 2 )) N [a'(ui, u 2 ) = 9 li^+t(i^Ti't)+i ] and D ' ~ (3/+1) 2 awa ^ from the 
critical point. We find a similar variation of J 7 (u, u + 5) [Fig. 7 (top) ] and D (u) [Fig. 7 (bottom) ] 
near the QCP u = 0, i.e., u = 00 as in the case of the QCP u = 0. T (u, u + S) falls sharply at u = 
and the fall becomes faster as we increase N. The quantity ln\D (u)\ increases at an enhanced rate 
and tends to blow up in the vicinity of u — as we increase the value of N. The inset of the figure 
shows that curves plotted in rescaled units collapse onto a single curve for different values of AT. 

The rescaled quantity D Jf is found to be a function of Nu 2 only. This feature of data collapse is 
analogous to the scaling behaviour of observables in the vicinity of a critical point. The finite size 
scaling hypothesis, valid in the critical region, is given by X^ = N~ Q{N \g — g c \ v ) where Xjy is 
some observable with the divergent behaviour Xm ~ \g — .9c|~ p close to the critical point g = g c . 
The exponent v is the correlation length exponent. In the present case p ~ and v = 2. 

IV. DISCUSSIONS 

In this paper, we have studied a S = \ spin ladder model the exact ground states of which are MP 
states. The ground state is spontaneously dimerized and doubly degenerate (broken translational 
symmetry) at all values of the parameter u excepting the points at u = and u = 00. At u = 0, a 
QPT occurs to the Haldane phase of an effective 5=1 chain which is signalled by the vanishing 
of a singlet excitation gap. The elementary singlet excitation in the dimerized phase is neither 
a magnon nor a spinon but a soliton in the dimer order. The lowest soliton excitations occur in 
pairs. At u = 00, there is another QPT to the rung-singlet phase with the vanishing of a triplet 
excitation gap, associated with triplet solitons. The ground states in all the three phases: Haldane 
(u = 0), spontaneously dimerized (0 < u < 00) and rung-singlet (u = 00) are spin liquids with 
no conventional long-range order in the two-spin correlation functions but are characterized by 
other types of order parameters. The spontaneously dimerized phase has long range order in dimer 
correlations which vanishes for u — > 0, 00 but there is no exponential tail. The Haldane phase 
has the string order parameter [HI Q21 H3 [14] whereas the rung-singlet phase has dimer-dimer 
correlations with the dimers located on the rungs. The two-spin correlation length is finite at u = 
and diverges as u — > 00 but no long range order develops in the latter case since the amplitude of 
spin correlations falls to zero in this limit. 

As pointed out in [16], QPTs in MP states are unconventional with the ground state energy 
analytic at g — g c , the transition point. A conventional QPT is signalled by a non-analyticity in 
the ground state energy. One can, however, generalize the definition of QPT to include cases where 
any observable quantity becomes non-analytic as the transition point is reached. MP states are 
an important class of states which provide an exact representation of many-body ground states of 
specific Hamiltonians. Also, every state of a finite system has an MP representation which thus 
provides the basis of the powerful DMRG method. In the thermodynamic limit, second order QPTs 
occur in MP ground states accompanied by vanishing energy gaps and diverging correlation lengths. 
We have studied the variation of the entanglement measures S(i) and S(i,j) as a function of u in 
the ground state of the spin ladder model with QCPs at u — and 00. The major goal of our study 
is to identify signatures of QPTs, if any, in the quantum information theoretic measures associated 
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with entanglement and fidelity. We provide a summary and analysis of our results below. 

Both S(i) and S(i,j) have zero values at u = oo, i.e., in the rung singlet phase (Figs. 2 and 
4) and nonzero values in the dimerized phase < u < oo. The entanglement measures can thus 
be treated as an order parameter with zero value at the QCP u — oo and non-zero value in the 
preceding dimerized phase. In the rung singlet phase, each rung is described by a spin singlet which 
is maximally entangled but the rung is disentangled from the rest of the system. The EL, £e, as 
calculated from S(i,j) diverges as u — > oo (Fig. 5 (bottom)) with £ B = ^, & being the spin-spin 
correlation length. The entanglement content in this case vanishes with infinite entanglement range. 
At the QCP u = 0, the entanglement measures have the magnitudes associated with the AKLT 
state of a spin-1 model. The entanglement measure has a local minimum at this point, rises to the 
maximum value at u — 1 and then decreases to the global minimum value zero at u — oo. The first 
derivatives of S(i) and S(i,j) both fall sharply to zero at u — and u — oo. The double derivatives 
of these quantities diverge as the QCPs are approached (insets of Figs. 1 — 4). The divergence 
arises from the structure of the von Neumann entropy involving terms such as log2U 2 or log^^s- 
A similar type of divergence occurs in the QPT of a model studied in [6]. We thus find that the 
entanglement measures S(i) and S(i,j) do develop distinctive features close to the QCPs u = 
and u = oo. 

We further looked for signatures of QPTs via the fidelity measure. Fidelity, i.e., the overlap of 
ground states for slightly different Hamiltonian parameters, is expected to drop abruptly at a QCP 
indicating a dramatic change in the ground state structure. We plotted F(u, S) = (u\u + S) with 
u and N for 6 = 10~ 3 and found that the quantity indeed falls to zero rapidly as the QCPs u = 
and oo are approached. The quantity ln\D(u)\, where D(u) is related to the second derivative of 
F, also provides a good signature of QPTs. You et al [32J has introduced a quantity, the so-called 
fidelity susceptibility \f which is defined as 

—2 In T(u, 5) 

Xf(u) = hms->Q p (25) 

One can easily check that \f has the same form as D(u). The finite size scaling hypothesis, which 
is expected to be valid in the vicinity of a QCP, leads to the collapse of curves onto a single scaling 
function (inset of Fig. 7) as the QCP u = oo is approached. The fidelity measures exhibit similar 
features in the case of a conventional QPT. The spin ladder model studied in the paper has spin 
liquid-type ground states with none of the phases exhibiting long range magnetic order in the two- 
spin correlation functions. The model has three distinct phases with characteristic quantum order 
parameters. A characterization of the transitions between the phases in terms of entanglement 
and fidelity measures provide a new perspective on the many body finitely correlated states and 
the transitions between them. Quantum information theoretic measures such as entanglement and 
fidelity provide a novel characterization of QPTs occuring in many-body condensed matter systems 
|33 | 134 } [35] . The present study illustrates this in the case of a spin ladder model with spin-liquid 
type ground states. 
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FIG. 2: Plot of S(i) (top) and (bottom) as functions of \. The inset (bottom) 

shows the diverging behavior of the second derivative of S(i) near u = oo. 
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FIG. 3: Plot of S(i,j) (top) and ^%p- (bottom) as functions of u for n = 1000. The 
inset (bottom) shows the diverging behavior of the second derivative of S(i,j) near 

u = 0. 
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FIG. 4: Plot of S(i, j) (top) and ^f^p (bottom) as functions of \ for n = 1000. The 
inset (bottom) shows the diverging behavior of the second derivative of S(i,j) near 

U = 00. 
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FIG. 7: Plot of ^(u, u + 6) (top) as a function of u and iV for 6 = .001 and in|D'(u)| 
(bottom) (N = 10 2 , 10 3 , 10 4 and 10 5 ) as a funct ion of u. The inset shows the data collapse 

for the rescaled function D Jf for same values of N . 
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